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Based on electronic structure calculations and theoretical analysis, we predict the (111) thin films 
of the SnTe class of three-dimensional (3D) topological crystalline insulators (TCIs) realize the 
quantum spin Hall phase in a wide range of thickness. The nontrivial topology originates from the 
inter-surface coupling of the topological surface states of TCIs in the 3D limit. The inter-surface 
coupling changes sign and gives rise to topological phase transitions as a function of film thickness. 
Furthermore, this coupling can be strongly affected by an external electric field, hence the quantum 
spin Hall phase can be effectively tuned under the experimentally accessible electric field. 

PACS numbers: 73.21.-b, 73.43.-f, 73.61.-r 


Topological crystalline insulators (TCIs) are topologi¬ 
cal phases of matter protected by crystal symmetries [1[ . 
TCIs have attracted wide attention after the prediction 
and observation of this phase in IV-VI semiconductors 
SnTe, Pbi_ x Sn x Se and Pbi_ x Sn x Te 2 [ 5 ] . The topologi¬ 
cal nature of the SnTe class of TCIs arises from the mir¬ 
ror symmetry of the rocksalt structure, and is manifested 
in the presence of topological surface states consisting 
of an even number of Dirac fermions. Recent experi¬ 
ments S0 have observed that two of the four Dirac 
fermions on the (001) surface become gapped when the 
structural distortion takes place and breaks one of the 
two mirror planes, as predicted by theory[2. [s[. This ob¬ 
servation establishes the topological protection by crystal 
symmetry—the defining property of TCIs. 

There is currently an intensive experimental investiga¬ 
tion of IV-VI TCIs in the low-dimensional and nanostruc¬ 
ture form[9l-[l6j. Recently we predicted that (001) films 
of TCIs with an odd number of atomic layers can pos¬ 
sess spin-filtered gapless edge states. Interestingly, the 
edges states can be gapped by a perpendicular electric 
field that breaks the mirror symmetry. Based on this 
novel functionality, we proposed a topological transistor 
device, in which spin and charge transport are coupled 
and can be simultaneously controlled by the electric field 
effect [l7| . 


Due to the cubic structure in the bulk, electronic struc¬ 
tures of IV-VI semiconductor thin films depend crucially 
on the growth direction [l8l-[2ll| . In this work, we study 
(111) thin films of IV-VI TCIs, which have been grown 
epitaxially in recent experiments [15, Q- We find that 
(111) thin films can host the quantum spin Hall state[22- 
l24j in a wide range of film thickness with both an even or 
odd number of layers (despite the different crystal sym¬ 
metries in the two cases). The largest nontrivial gap is 
found to be around 50 meV. Here the quantum spin Hall 
phase arises from the coherent coupling between the top 
and bottom topological surface states of TCIs in the 3D 
limit. This inter-surface coupling mechanism has been 
previously proposed in model studies of B^Ses thin films 


25-i3|, which inspired part of this work. Furthermore, 
we find that the inter-surface coupling in (111) TCI thin 
films can be strongly tuned by a moderate external elec¬ 
tric field, leading to an electrical control of the topological 
property. 

For convenience, we will only show below the calcu¬ 
lations for SnTe as a representative of TCIs in IV-VI 
semiconductors. All the qualitative conclusions apply 
to Pbi_ x Sn x Te and Pbi_ x Sn x Se as well. The rock- 
salt structure of SnTe has two inequivalent (111) sur¬ 
faces, ending at the plane of Sn or Te atoms respec¬ 
tively. Both surface states consist of four Dirac cones 
centered at time-reversal-invariant momenta: one is at 
T point and the other three are at M points [ 28 U 3 H . 
as observed in recent ang le-resolved photoemission spec¬ 
troscopy experiments [l4 32 . |j3]. Different from (001) 


surface states, these four Dirac cones here are not all 
equivalent. Specifically, the Dirac point at T point is well 
inside the bulk gap while the Dirac points at M points 
are very close to the valence or conduction band edges. 
Later, we will show this feature has important conse¬ 
quences for the electronic structures and band topology 
of (111) thin films. 

In a thin film, the top and bottom surface states hy¬ 
bridize to open up energy gaps at these four Dirac points. 
It is important to treat films with an odd and even num¬ 
ber of layers separately due to the difference in crys¬ 
tal symmetry: the former has inversion symmetry and 
the latter doesn’t. Our calculations based on the tight- 
binding model of Lent et al [HU show different band struc¬ 
tures for the two cases, as shown in Fig.l. Nonetheless, 
both are semiconductors with a small fundamental gap 
located at or near the T point, and a much larger band 
gap at M. This is because the penetration length of sur¬ 
face states at M is much larger than the one at T, hence 
the hybridization-induced gap at M is correspondingly 
larger than the one at T by orders of magnitude. As a 
result of this inequivalence between T and M, (111) films 
of TCIs in the SnTe class are narrow-gap semiconductors, 
whose properties are controlled by low-energy physics at 
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FIG. 1: The typical band structures of (111) thin films: 24- 
layer (top) and 25-layer (bottom). The different symmetries 
of even and odd number of layers thin films (see insets) leads 
to their qualitatively different band structures. 

the T point. 

To deduce the band topology, we build an effective k-p 
Hamiltonian at T. First considering the case of an even 
number of layers, here the top and bottom surface are 
different, terminated by Sn-atom and Te-atom respec¬ 
tively, which breaks the inversion symmetry (see insets 
in Fig. 1). Each surface is described by a two-dimensional 
Dirac fermion: H t ,b = ^t,b(k x s y — k y s x ) + E t ,b, where v t 
and Vb denote the Dirac velocities of the top and bottom 
surface respectively; E t and Eb denote the Dirac point 
energies of the two surfaces. Denoting the two surfaces 
by r z = ±1, the hybridization term to zeroth order in k 
is given by r x . This leads to a four-band k-p Hamiltonian 
for the film: 

-H(k) — v(k x s y k y s x )T z T T (1) 

Here for simplicity we have assumed that v t = — Vb = v: 
including two velocities of different magnitude will not 
change the conclusion of band topology determined from 
(PQ). S = ( E t — Eb)/2 measures the energy difference of 
the Dirac points at top and bottom surface, which can 
be affected by a perpendicular electric field. 

By diagonalizing the k), we obtain four bands with 
energy-momentum dispersion Ey( k), Ew{ k), —Ey(k) 
and —Ew{ k), where £V,w(k) are given by 

Ev,w{ k) = \Jm 2 + (yibj; + k%)v 2 ±"<5)2 (2) 

The (5-term causes the two degenerate Dirac cones split 
in energy. For m = 0, the hole-band of higher Dirac cone 
and the electron-band of lower Dirac cone will cross each 
other at E = 0, forming an circular Fermi surface defined 
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FIG. 2: The effect of thickness (top) and the effect of ex¬ 
ternal perpendicular electric field on 12-layer thin him (bot¬ 
tom). The quantum spin Hall (QSH) phase can be achieved 
by changing the thickness or tuning the external electric held. 
The red/blue dots in the inset of top panel represent that the 
states mainly compose of Te or Sn atoms and the size means 
the corresponding weight. 

by (k 2 + ky)v 2 = S 2 . Now turning on the hybridization 
term m/0 will open a gap at this circular Fermi surface: 
Eg = 2 1 ra |. In addition, the gap at T point increases to 
E g ( r) = e l\Jm 2 + S 2 . The resulting band structure is 
particle-hole-symmetric and has Mexican-hat dispersion, 
with degenerate band gaps around a circle in k space. 
This is indeed found in our band structure calculation 
shown in the top panel of Fig. 1. 

Importantly, the band structure m shows that the 
band gap does not close when S is tuned to zero, as long 
as m stays hnite. Therefore, the Z 2 topology can be 
determined by setting (5 = 0, and is controlled entirely 
by m, the coupling between top and bottom surfaces. 
The Hamiltonian for S = 0 takes the form of a two- 
dimensional massive Dirac Hamiltonian, as in Bernvig- 
Hughes-Zhang model for HgTe/CdTe quantum well[23|. 
The sign of m controls the Z 2 topology. Similar to previ¬ 
ous proposals for Bi 2 Se 3 thin films [25H271]. we expect that 
the inter-surface hybridization m, both its magnitude 
and sign, should depend on the film thickness, thereby 
leading to gap change and topological phase transition 
as a function of film thickness. 

This is confirmed by our tight-binding calculations. As 
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FIG. 3: The edge states of 12-layers thin film (left panel), 14- 
layer thin film(middle panel) and 12-layer thin films under 
electric field (right panel). For 14-layer thin film, there are 
gapless states with linear dispersion exist in the bulk gap, 
which confirm that 14-layer thin film is in the quantum spin 
Hall (QSH) phase. Such kind of edge states are absent for 12- 
layer thin film, a normal insulator (NI), while we can drive 12- 
layer thin film into QSH phase by external electric field. (The 
electric potential between Te-terminated and Sn terminated 
surfaces is -0.2 V) 


shown in top panel of Fig. 2, both the fundamental gap 
Eg and the gap at Y point E g (Y) change with the film 
thickness. When the thickness deceases, both gaps ex¬ 
hibit a oscillatory behavior, confirming our expectation 
from the above k • p model. The gap oscillation is ex¬ 
pected to signal sign change of the Dirac mass m, lead¬ 
ing to Z 2 topological phase transitions between 12-layer 
and 14-layer thin films. Furthermore, as shown in the 
inset of top panel of Fig. 2, the conduction bands of 12- 
layer thin film mainly compose of Sn orbitals, while it’s 
mainly from Te orbitals for 14-layer thin film, suggesting 
that 12-layer and 14-layer thin films belong to different 
topological phases. 

To determine which thickness range is topological non¬ 
trivial, we directly calculate the edge states in an semi¬ 
infinite geometry by using the recursive Green’s func¬ 
tion method [35]. Without loss of generality, we take 
armchair-type edges as an example. As shown in Fig. 
3, for 14-layer thin film, the linear-dispersion gapless 
states exist in the bulk gap with a two-fold degeneracy 
at T point, named Kramer’s degeneracy, forming a one¬ 
dimensional Dirac point. In comparison, there are no 
such kind of helical edge states in the 12-layer thin film. 
From such a calculation for other film thickness, we estab¬ 
lish unambiguously a topological phase diagram shown in 
Fig. 2. In particular, all the even number of layers thin 
films between 14 and 30 layers are quantum spin Hall 
insulators. 

In addition to its sensitivity on the thickness, we find 
that the inter-surface coupling m can also be strongly 
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FIG. 4: The phase diagram for the even number layer of thin 
films varying with thickness and external electric field. 

affected by a perpendicular electric field [36[. Here the ef¬ 
fect of an electric field is modeled by a linearly-increasing 
electrostatic potential across the film thickness direction, 
leading to a potential difference (j) between the Te- and 
Sn-terminated surface. We carefully calculate the evolu¬ 
tion of the band structure of 12-layer thin film by vary¬ 
ing (j) and find the gap exhibits an very interesting os¬ 
cillated behavior: it indeed closes at some critical val¬ 
ues, implying a topological phase transition induced by 
the electric field. Specifically, when we apply a nega¬ 
tive electric field, the gap will first decrease and closes at 
(j) = —0.079V. As the electric field further increases, the 
gap reopens and achieve the maximum around 75 meV 
around 0 — —0.51V and then goes down. If the electric 
field is strong enough, the thin film will become to be a 
metal. It is worth noting that the critical electric field 
strength for topological phase transition is only around 
4 x 10 7 V/m, which is an order smaller than the typical 
breakdown field strength for SnTe related thin films(37j. 
To further verify this explicitly, we calculated the edge 
states for 12-layer thin film under different electric field. 
As shown in the Fig. 3, without electric field, the edge 
states are fully gapped for the 12-layer thin film, while 
when the electric field is stronger than <f> = —0.079V, 
there will be topologically protected edge states in the 
bulk gap (a typical results with </> = —0.2V shown in the 
right panel). Similar behavior is seen at positive electric 
fields, with topological phase transitions at (j) = 1.1V and 
2.2V. The corresponding phase diagram is shown in Fig. 
2 . 

By doing similar calculations for other film thickness, 
we obtain the complete phase diagram as a function of 
thickness and electric field, as shown in Fig. 4. From this 
we find that the quantum spin Hall phase of thin films 
thicker than 16 layer is very robust, whereas the thin films 
between 10-layer and 14-layer are very close to the phase 
boundary. Topological phase transition can be realized 
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FIG. 5: The subband energy level diagram at F point and the 
topological phase diagram for odd number of layers (111) thin 
films with Sn-atom (top panel) and Te-atom (bottom panel) 
terminated. The energy order of T + and T_ states changes in 
oscillatory behavior with the number of layers, which results 
in corresponding the topological phase oscillation. The Dirac 
mass m (gap at T point) is m = F?(r+) — (T_). 


under an experimentally feasible external electric field, 
which can be set up by using top or bottom gate (36). 

Now we turn to the case of the odd number of layers 
thin films, where both surfaces are terminated by the 
same atoms (Sn or Te atoms) and inversion symmetry is 
present (see the insets in Fig. 1). The energy difference 
term Sr z in the k-p Hamiltonian (pQ) is no longer allowed, 
leading to the standard Dirac Hamiltonian H = ma x + 
v(k x Sy — k y s x )a z . The presence of inversion symmetry 
allows us to determine the topological invariant Z2 based 
on the parity cr iter ion [38j. Clearly, the conduction and 
valence bands at T are respectively the bonding and anti¬ 
bonding hybridization of top and bottom surface states, 
with opposite parities. Tight-binging calculations show 
that the 2D band structure near F is described very well 
by a single-valley Dirac fermion with a small mass, and 
the dispersion is linear and particle-hole symmetric in a 
sizable energy range (Fig.l). 

Fig.5 shows the energies of subbands at T for films 
thinner than 95 layers. The Dirac mass m, given by 
the energy difference between the bottom conduction and 
top valence band labelled by T + and T_, depends non- 
monotonously on the thickness. Importantly, m changes 
sign at certain critical thickness: 19-, 39-, 63- and 85- 


layer for Sn termination, and 13-, 27-, 47, 63- and 83-layer 
for Te termination. The Dirac mass reversal switches the 
parity of conduction and valence band edges, resulting 
in topological phase transition between a trivial insula¬ 
tor and a quantum spin Hall insulator. According to our 
tight-binding calculation, the band gap in the quantum 
spin Hall phase can reach 18meV at 25 layers for Sn ter¬ 
mination, and 66meV at 19 layers for Te termination. 
The topological phase transition tuned by film thickness 
is similar to HgTe/CdTe quantum well[23|. Similar to 
even number of layers thin films, perpendicular electric 
potential can also induce topological phase transition for 
the odd number of layers thin films. 


The mechanism for realizing quantum spin Hall phase 
in TCI thin films via inter-surface coupling is similar to 
the thin films of topological insulators B^Ses Hi-ill, 
which partly motivated this work. Due to the much 
larger penetration length of SnTe surface states, the 
hybridization gap is much larger than that in B^Ses. 
Furthermore, the proposed quantum spin Hall phase in 
Bi 2 Se 3 was found to depend very sensitively on the film 
thickness (251427) . with Z 2 topology changing back and 
forth with a period of only one or two quintuple layers. 
In comparison, our results show that the quantum spin 
Hall states in TCI films exists in an extended thickness 
range and hence are much more robust. Finally, our theo¬ 
retical analysis of electronic and topological properties of 
TCI thin films based on inter-surface hybridization may 
apply to many other TCI materials that are currently 
being studied (391-143). 

Note added: After the completion of this work[44(, we 
learned of an independent work on (111) thin films of 
TCIs without the external electric field [45|. 


We thank Tim Hsieh for helpful discussions, and R. 
Buczko for alerting us to his work. This work was sup¬ 
ported by the STC Center for Integrated Quantum Ma¬ 
terials, NSF Grant No. DMR-1231319. 


[1] L. Fu, Phys. Rev. Lett. 106, 106802 (2011). 

[2] T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil and L. 
Fu, Nat. Commun. 3, 982 (2012). 

[3] Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, 
T. Takahashi, K. Segawa, Y. Ando, Nat. Phys. 8, 800 
( 2012 ). 

[4] P. Dziawa, B. J. Kowalski, K. Dybko, R. Buczko et al ., 
Nat. Mater. 11 , 1023 (2012). 

[5] S.-Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian et 
a/., Nat. Commun. 3, 1192 (2012). 

[6] Y. Okada, M. Serbyn, H. Lin, D. Walkup et al. , Science, 
341, 1496 (2013). 

[7] I. Zeljkovic, Y. Okada, M. Serbyn, R. Sankar, D. Walkup, 
W. Zhou, J. Liu, G. Chang, Y. J. Wang, M. Z. Hasan, 
F. Chou, H. Lin, A. Bansil, L. Fu, V. Madhavan, 
arXiv: 1403.4906 (2014) 

[8] M. Serbyn and L. Fu, Phys. Rev. B 90, 035402 (2014). 








































5 


[9] M. Safdar, Q. Wang , M. Mirza , Z. Wang , K. Xu , and 
J. He, Nano Lett. 13, 5344 (2013) 

[10] A. A. Taskin, F. Yang, S. Sasaki, K. Segawa, and Y. 
Ando, Phys. Rev. B 89, 121302 (2014) 

[11] J. Shen, J. J. Cha, Nanoscale 6, 14133 (2014) 

[12] J. Shen, Y. Jung, A. S. Disa, F. J. Walker, C. H. Ahn, J. 
J. Cha, Nano Lett. 14, 4183 (2014) 

[13] H. Guo et al, APL Materials 2, 056106 (2014) 

[14] C. Yan et al. , Phys. Rev. Lett. 112, 186801 (2014) 

[15] B. A. Assaf, F. Katmis, P. Wei, B. Satpati, Z. Zhang, S. 
P. Bennett, V. G. Harris, J. S. Moodera and D. Heiman, 
Appl. Phys. Lett. 105, 102108 (2014) 

[16] S. Sasaki, Y. Ando, arXiv:1410.4852 (2014) 

[17] J. Liu, T. H. Hsieh, P. Wei, W. Duan, J. Moodera, L. 
Fu, Nat. Mat. 13, 178 (2014) 

[18] R. Buczko and L. Cywinski, Phys. Rev. B 85, 205319 

( 2012 ). 

[19] X. Qian, L. Fu, J. Li, arXiv: 1403.3952 (2014) 

[20] H. Ozawa, A. Yamakage, M. Sato, and Y. Tanaka, Phys. 
Rev. B 90, 045309 (2014) 

[21] M. Ezawa, New J. Phys. 16 065015 (2014) 

[22] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005); ibid , 95, 146802 (2005). 

[23] B. A. Bernevig, T. L. Hughes and S. C. Zhang, Science 
314,1757 (2006). 

[24] M. Konig, S. Wiedmannl, C. Briine, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X.-L. Qi, S.-C. Zhang, Science 
318, 766 (2007). 

[25] J.Linder, T. Yokoyama, and A. Sudb0, Phys. Rev. B 80, 
205401 (2009). 

[26] C. X. Liu, H. J. Zhang, B. Yan, X. L. Qi, T. Frauenheim, 
X. Dai, Z. Fang, and S. C. Zhang, Phys. Rev. B 81, 
041307 (2010). 

[27] H.-Z. Lu, W.-Y. Shan, W. Yao, Q. Niu, and S.-Q. Shen, 
Phys. Rev. B 81, 115407 (2010) 

[28] J. Liu, W. duan and L. Fu, Phys. Rev. B 88, 
241303(R)(2013) 

[29] J. Wang, J. Liu, Y. Xu, J. Wu, B. L. Gu, W. Duan, Phys. 
Rev. B 89, 125308 (2014) 

[30] Y. J. Wang, W.-F. Tsai, H.Lin, S.-Y. Xu, M. Neupane, 
M. Z. Hasan, and A. Bansil, Phys. Rev. B 87, 235317 
(2013) 

[31] S. Safaei, P. Kacman, and R. Buczko, Phys. Rev. B 88, 
045305 (2013) 

[32] Y. Tanaka, T. Shoman, K. Nakayama, S. Souma, T. Sato, 
T. Takahashi, M. Novak, K. Segawa, and Y. Ando, Phys. 
Rev. B 88, 235126 (2013) 

[33] C. M. Polley et al. , Phys. Rev. B 89, 075317 (2014) 

[34] C. S. Lent, M. A. Bowen, J. D. Dow, R. S. Allgaier, O. F. 
Sankey and E. S. Ho, Superlattices Microstruct. 2, 491 
(1986). 

[35] M. P. L. Sancho, J. M. L. Sancho, J. M. L. Sancho and 
J. Rubio, J. Phys. F: Met. Phys. 15, 851 (1985) 

[36] X. Qian, J. Liu, L. Fu and J. Li, Science 346, 1344 (2014) 

[37] S. Parvanov, V. Vassileva and L. Aljihmania, J. Opto- 
electron. Adv. Mater. 7, 1299 (2005) 

[38] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 
39 M. Kargarian and G.A. Fiete, Phys. Rev. Lett. 110, 

156403 (2013). 

[40] M. Kindermann, arXiv: 1309.1667 (2013). 

[41] M. Ye, J. W. Allen, K. Sun, arXiv:1307.7191 (2013) 

42 H. Weng, J. Zhao, Z. Wang, Z. Fang, and X. Dai, Phys. 

Rev. Lett. 112, 016403 (2014). 

[43] T. H. Hsieh, J. Liu, and L. Fu, Phys. Rev. B 90, 


081112(2014) 

[44] J. Liu, et al ., http://meetings.aps.org/Meeting/MAR13/ 
Event/185363 (2013) 

[45] S. Safaei, M. Galicka, P. Kacman and R. Buczko, 
arXiv: 1501.04728 (2015) 


